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We derive analytic solution for pulsed frequency conversion based on electromagnetically induced 
transparency (EIT) or maximum coherence in resonant atomic vapors. In particular drive-field 
and coherence depletion are taken into account. The solutions are obtained with the help of an 
Hamiltonian approach which in the adiabatic limit allows to reduce the full set of Maxwell-Bloch 
equations to simple canonical equations of Hamiltonian mechanics for the field variables. Adiabatic 
integrals of motion can be obtained and general expressions for the spatio-temporal evolution of 
field intensities derived. Optimum conditions for maximum conversion efficiency are identified and 
the physical mechanism of nonlinear conversion in the limit of drive-field and coherence depletion 
discussed. 



I. INTRODUCTION 

Resonant nonhnear optics in atomic gases has received 
a great impetus in recent years due to new concepts 
based on the apphcation of specific coherence and in- 
terference effects. Two of these concepts attracted par- 
ticular attention. The first mechanism uses the effect 
of electromagnetically induced transparency (EIT) |j, 
the second maximum coherence between two metastable 
atomic levels [Q. Both schemes allow for high-efficiency 
nonlinear conversion with substantially alleviated phase- 
matching problems in a rather dilute ensemble of atoms 
if the medium is driven by an undepleted coherent cou- 
pling field or if an undepleted atomic coherence is as- 
sumed. These assumptions do not take into account the 
resources needed to maintain the drive field or atomic 
coherence. Thus in the considered limit the overall 
efficiency of the resonant nonlinear processes is small. 
In the present paper we analyze and compare EIT- 
and maximum-coherence based systems with each other 
and with the conventional non-resonant schemes of non- 
linear optics taking into account the transparency- or 
coherence-maintaining fields and their depletion. Effects 
of inhomogeneous broadening are however disregarded. 
Deriving analytic solutions for the nonlinear pulse prop- 
agation problem, we show that in both systems also a 
large overall conversion efficiency can be achieved. The 
physical nature of the conversion can however no longer 
be associated with EIT or maximum coherence alone. 

In EIT a strong resonant electromagnetic (em) field 
of frequency uj2 is applied to the transition between two 
excited states |2) and |3), (Fig. la), causing a splitting of 
both (Autler-Townes effect). When a weak probe field 
of frequency L03 resonant with the transition from the 
ground-state |1) to the bare state |3) is applied, its lin- 
ear interaction with the medium is almost perfectly can- 
celled. The vanishing of the linear susceptibility, i.e. ab- 
sorption and refraction, for the probe field is due to a 
destructive interference between the two excitation paths 
through the Autler-Townes doublet and persists even for 
small splittings. If the ground state |1) is coupled to the 
intermediate state |2), on the other hand (e.g. by a two- 
photon transition), photons of frequency ^3 are gener- 



ated. While the resonant linear absorption of these pho- 
tons is suppressed by destructive quantum interference 
(EIT), the nonlinear susceptibility responsible for the 
generation of them is only affected by the splitting, which 
can be rather small. Since there are no resonant contribu- 
tions to the refractive index, there is furthermore perfect 
phase matching except for contributions from other off- 
resonant transitions. Owing to the large dispersion near 
the transparency frequency of EIT, a small detuning 62 
can be applied to compensate for these contributions and 
phase matching can easily be obtained without sacrific- 
ing the cancellation of absorption. Thus EIT provides a 
perfect system for large nonlinear conversion with a min- 
imum of atoms. Since the first demonstration of phase 
matching in EIT-assisted four- wave frequency mixing , 
many other experiments confirmed considerable improve- 
ment in the conversion efficiency when EIT has been used 
(see, e.g., Refs. [§ |). 

Another mechanism proposed recently is referred to 
as nonlinear optics with maximum coherence |^. The 
idea is here to prepare and maintain the atoms in a coher- 
ent superposition of atomic states |1) and |2) with equal 
amplitudes. This can be done either by two strong fields 
exciting the |1) — |2) coherence via a Raman transition 
or by other means such as rapid adiabatic passage [D . 
The coherence, established in such a way, plays the role 
of a strong local oscillator. If then a relatively weak field 
L02 detuned by a sufficiently large amount S3 from the 
|2) — |3) transition is applied. Fig. lb, it will beat against 
the local oscillator to generate the sum (or difference) 
frequency. Here, the refraction due to resonant transi- 
tions does not vanish. But since the nonlinear coupling 
strength is of the same order as the linear susceptibility, 
complete conversion occurs over a distance smaller than 
the coherence length Consequently there is no need 
to phasematch the propagating beams. We note that the 
relation between the |1) — |2) coherence and the wave 
mixing processes has been pointed out in earlier work 
on nonlinear optics (see, e.g., Refs. |^ and references 
therein) . 

In the limit of undepleted drive or constant coherence, 
the nonlinear conversion process in both the EIT and 
maximum coherence scheme affects the quantum state of 
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the atoms only in the perturbative sense. Consequently 
the probability of photon conversion per atom is small. 
For the same reason the overall conversion efficiency in- 
cluding the transparency and coherence generating fields 
is small. We here analyze the properties of the two mech- 
anisms for arbitrary strength of the drive fields taking 
into account their depletion and compare them to con- 
ventional off-resonant nonlinear optics. 

The drive depiction turns the propagation problem 
into a truly nonlinear one. Its solution is particularly 
challenging for pulses and is in general possible only nu- 
merically. In order to obtain transparent analytical solu- 
tions we here apply the so-called Hamiltonian approach 
[T^ , p^ , [T^ which allows for a solution in a wide range 
of physically relevant situations. This approach is espe- 
cially useful under adiabatic conditions, i.e. when the 
atoms are excited by the laser pulses in such a way that 
they remain in the same instantaneous eigenstate of the 
interaction Hamiltonian during the entire process. In the 
present work, we will assume that the interaction is adi- 
abatic. The adiabatic approximation requires a slow rate 
of evolution as compared to the frequency separation of 
the adiabatic eigenstates. This results usually in a re- 
quirement for the product of the pulse duration and the 
Rabi frequency of the radiation field to be much larger 
than unity. Thus, for sufficiently intense fields, the pro- 
cess can be adiabatic even for short pulses. It should be 
noted however that this assumption rules out such im- 
portant effects as group velocity reduction which results 
from lowest-order non-adiabatic corrections. For suffi- 
ciently large intensities or cw fields these effects do not 
influence the process. An analysis of non-adiabatic cor- 
rections to the Hamiltonian approach taking into account 
group delays will be given in a future publication. 

To simplify the calculations and the interpretation of 
the analytic results, we restrict ourselves to a somewhat 
idealized three-level atomic excitation scheme (Fig. 1). 
The nonlinear three-wave mixing in a three-level system 
is normally forbidden due to symmetry However, this 
scheme provides the simplest example where all above 
mentioned mechanisms of nonlinear-optical wave mixing 
may take place. It should be noted that similar results, 
however with much lengthier expressions and larger num- 
ber of atomic parameters, can be found for a more re- 
alistic three-level scheme where |1) — |2) transition is a 
two-photon one psf . Moreover, the 3- wave-mixing pro- 
cesses is possible when a dc electric field is applied to the 
atomic sample p6| . 

The paper is organized as follows. In Sec. II we dis- 
cuss resonant nonlinear optical processes based on EIT 
or maximum coherence for undepleted drive fields or un- 
depleted coherence, respectively, and compare them to 
conventional off-resonance nonlinear optics. In Sec. HI we 
outline the Hamiltonian approach, which allows to elim- 
inate the atomic degrees of freedom assuming adiabatic 
following and to map the pulse propagation problem to 
the dynamics of a one-dimensional nonlinear pendulum. 
Making use of this formalism we derive full analytic so- 



lution for EIT and maximum-coherence based nonlinear 
optics in Sec. IV, taking into account drive-field and co- 
herence depletion. 

II. RESONANT NONLINEAR OPTICS IN A 
COHERENTLY PREPARED 3-LEVEL SYSTEM 

We first consider the propagation of pulsed em fields 
in a medium of three-level atoms (Fig. 1) in which either 
a constant drive field mixes the two excited states |2) and 
1 3) or in which a constant coherence between the lower 
two states |1) and |2) is maintained. The ffist case corre- 
sponds to resonant nonlinear frequency conversion based 
on EIT D, |1, the second one to nonlinear optics with 
maximum coherence The electric field propagating 
in the z-direction is assumed to consist of three compo- 
nents with carrier frequencies wi, 0^2 and L03 = loi -\- lo2'. 

E{z, t) — ^|^£j(z, t) exp(— i(tjji — kjz)) -I- c.c^j . (1) 
j 

Here kj = njuj^jc with rij being the refractive index at 
frequency LOj due to levels outside the three-level system 
of Fig. 1. This background refraction gives rise to the 
"residual" phase mismatch determined as: 

Afc = /ci + ^2 - k-i. (2) 
(a) (b) 



I3> |3> 




FIG. 1: Resonant sum-frequency generation, (a) strong drive 
field (LJ2) between metastable state |2) and excited state |3) 
creates EIT; (b) strong drive lo\ (two-photon, Raman or mag- 
netic coupling) generates maximum coherence between |1) 
and 1 2). 

In the approximation of slowly varying amplitudes and 
phases Maxwell's propagation equations read in a moving 
frame 



where Ej and Vj are functions of the coordinate z and 
the retarded time t = t — z / c. Vj are the components of 
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the medium polarization: 

P = "^^(j^j cxp(— — kjz)) + c.c.^ , 
j 

which can be expressed in terms of the atomic probabihty 
amphtudes c„ in levels |1), |2) and |3): 



Vi = Ndiclc2, 



(4) 
(5) 
(6) 



N being the density of active atoms, di, d2 and d^ are 
the real dipole moments of the transitions |1) |2), 
|2) |3) and |1) |3) respectively. 9 is phase of the 
|2) — > |3) dipole moment, which in general cannot be 
chosen freely. 

Assuming decay only from the top most state |3) with 
rate 7 the state amplitudes obey in rotating wave ap- 
proximation the equations 



ci = ifli C2 + i^'^ C3, 

C2 = inici+ini.e''^'"^''^ C3 + iS2C2, 



(7) 



C3 



in2 ^-'"+'^'''-02 + i^sCl + I (53 + ll) C3, 



where ^2 and ^3 are the frequency detunings indicated in 
Fig. 1. 



(52 = — W21 



(8) 



with LDij denoting the transition frequencies between the 
corresponding levels, fii, VL2 and ^3 are the Rabi fre- 
quencies for transitions |1) — |2), |2) ~ |3) and |1) — |3), 
respectively: 



2h 



(9) 



A. EIT with undepleted coupling field 

Here, we consider three fields interacting with the 3- 
level system. For simplicity we assume ^3 = 0. Further- 
more we consider the case of a strong, undepleted drive 
field with frequency UJ2 

\n2\ » |r!i|,|r!3l,7,|52|. 

Then, the solution of the atomic equations of motion (Q) 
with the initial condition ci {t — > —00) = 1 yields 

|ci|«l. 

Assuming quasi-adiabatic evolution, i.e. not too fast 
changing fields, we find 



C2 



C3 = 



\^2\'^~iS2l 
^1^26^''^+"^''' - ^362 



162! 



Substitution into Maxwell's propagation equations gives: 
9f 1 'KNd\ wi 7 



dz 



fl C 1^2 



n\n2\ c 

863 . TrNdl UJ3 62 



£3 



dz 



72 '^3 



fi C |f72| 

TTNdids UJ3 te+iAk: 

n\n2\ c 



(10) 



(11) 



Si. 



Equations (^^ and (^_1|) are linear differential equations, 
which can easily be solved. We now consider the case in 
which no £3 field is incident on the medium, £3 (^ = 0) = 
0. Introducing the normalized intensity (photon flux) 



Vj = 



I3 _ c\£,\ 



and the coupling strength 



2TTUJjd'j 



he 

the solution of equations (^0|) and ( [III ) reads: 



(12) 



(13) 



^ 2 e"^"' sin^ ( Kz-y/1 + {Ak'/2n) 



1 + {Ak'/2K 



(14) 



where 7710 = 771 (z = 0) is the photon flux at the entrance 
to the medium and we have introduced the conversion 
coefficient k: 



2 |r!2| ' 



the loss coefficient F: 



and 



7 

1^21 



Afc' = Ak 



Ms' 
N ^1362 

2 m^f 



(15) 



(16) 



(17) 



is the total phase mismatch, including the background 
value AA; and the contribution from resonant transition 
|1) |3). In the EIT case, the resonant contribution to 
the phase mismatch is always smaller than the conversion 
coefficient k by a factor S2/\^2\- Moreover, a small 
detuning 82 can be introduced to compensate the residual 
phase mismatch. Thus, the EIT scheme represents an 
ideal situation with complete phase matching, where the 
optimum conversion occurs for a density-length product 



Nz 



I opt 



(18) 
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In order to minimize the absorption losses, Eq. (16), 
the Rabi frequency of the couphng field |Jl2| has to be 
sufficiently large: 




For atomic media which we consider here, the off- 
resonant (background) contributions to the refractive in- 
dex Uj are of the order of j|] : 



(25) 



B. Undepleted coherence pi2 

In the following we discuss the scheme of resonant non- 
linear optics with maximum coherence. For this situa- 
tion, constant state amplitudes c\ and C2 are assumed, 
maintained e.g. by a constant strong drive field fii. In 
this case the wave vector of the atomic coherence c*C2 
is equal to k\ which can be changed by, e.g., the small 
detuning 82 |Q . The backaction of the atoms to this field 
is disregarded and thus the corresponding coupling does 
not need to be taken into account. Under the condition 

1(53 + «7l » \^i\ 

the amplitude of the excited state can be adiabatically 
eliminated, which yields 



C3 



^3 -I- 17 



(19) 



Substitution of Eqs. (^, (^, (19) into the Maxwell prop- 
agation equations (||) gives: 



9£2 
dz 



dz 



* h c (5| 72 



«7 I a r. 

C2 C2 



c O3 -f 7-^ 



iAkz 



nNdj CJ3 63 



c (5| -f 72 



«7 , ,2 r. 
Cl\ £3 



(20) 



(21) 



TTNd2d3 UJ3 S3 — i"f 
r — oPi2e 



iO+iAkz c 
2^, ^Pik:'- ^2: 



'% C 5|-f7 

where pi2 = \clc2\- 

This is again a set of linear differential equations whose 
solution has the same form as that for the KIT case, Eq. 
(p^), with the substitution rjiQ 7720 and corresponding 
parameters (assuming ^3 ^ 7): 



P12, 



r = 



7 P2 \C2\^ + M3 |ci|^ 



<53 

Afc' = Afc 



Pi2y/mla. 

iV^3|ci| ~ P2\C2\ 



(22) 
(23) 
(24) 



The total phase mismatch A/c' includes the background 
value Afc and the contributions from resonant transitions 
|1) ^ |3) and |2) -> |3). 



where pjm and 5jm are the coupling constants and de- 
tunings, respectively, for the wave with frequency ujj and 
transition to the far-detuned state \m) not belonging to 
the three-level system of Fig. 1. Since Sj^ ^ ^3, the 
resonant contributions to the phase mismatch are the 
dominant ones. As can be seen from Eqs.(p2|) and (|24|), 
they are in general of the same order as the conversion 
coefficient n if the atomic coherence pi2 is large (of the 
order of 1/2). Therefore, efficient energy transfer from 
the UJ2 field into the 1^3 field: 773(0 ~ 7720, occurs already 
within a length I = ~ Lc = 2/Afc', the coherence 
length. This feature constitutes the main advantage of 
the maximum coherence scheme over conventional non- 
linear optics, because the phase matching is no longer 
important. 

If the conversion length I is much smaller than the 
coherence length Lc, l/Lc <C 1, and the losses are 
small, (r/Z)2:|opt ^ 1, maximum conversion occurs for 
a density-length product 



Nz\ ^^Im^n^ . 



(26) 



l/Lc^ 1 is realized when the parameters are chosen such 
that Afc' is small, or even better, vanishes, i.e. if 

2Afc M2 |c2|^ - Ms |ci|^ 



N 



This can be achieved in different ways: (i) by tuning the 
wave vector fci of the atomic coherence (e.g., by intro- 
ducing a detuning 82, as in Ref. |^), (ii) by selecting 
the appropriate detuning ^3, as in Ref. 0, and/or, by 
preparation of atoms in a superposition with suitable am- 
plitudes Ci , C2 . 

In order for the absorption losses to be negligible 
within the optimum propagation distance, the following 
condition has to be satisfied: 

7 1 

03 Pl2 

which indicates once again that it is indeed advantageous 
to prepare a large atomic coherence pi2. 



C. Conventional nonlinear optics: weak excitation 

Now, we discuss the case of a weak excitation of atoms 
which corresponds to the regime of conventional nonlin- 
ear optics. The weak excitation takes place, for example, 
when both detunings are very large: 
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|J3|,|52|»|f^,|,7- 

Under this condition, the atomic probabihty amphtudes 
are: 



|ci| 

C2 

C3 



1, 



52 (^3 + il) 



Since |c2|,|c3| |ci|, the medium polarization at fre- 
quency L02 is much smaller than for the cji and 0^3 com- 
ponents: 17^21 ^ I'^'ihl'^al- Therefore, we can assume 
that the 0^2 field is almost undepleted, \^2 \ ~ const{z). 
The propagation equations read for this case: 



d£i 

dz 



d£3 

dz 



+i 



'^1 

fi-02 C 

nNdid^LUi 15^2 1 
he 6283 
. irNdl UJ3 S3 - ij 
h c i5| -I- 7^ 
.TrNdidsuJs \n2 



(27) 



AO — iAkz 



£3, 



£3 



(28) 



he 



52^3 



-ie+iAkz 



£i- 



The solution of this set of propagation equations has ex- 
actly the same form as in the maximum-coherence and 
EIT cases, Eq. (p^, with parameters: 



N 
2 


^2\^/^J.lfJ'3 
S2S3 


(29) 


7 

^3^ 


[Jh 52 

\l ^J'3 \^2y 


(30) 


Afc 


2 ^3 52)- 


(31) 



r 

2 \ 63 62 , 

We see that the conversion length I — is here much 
larger than both in the maximum coherence, Eq. (22), 
and in the EIT, Eq. (p^), schemes. In fact, the conversion 
length is also much larger than the (non-phase-matched) 
coherence length Lc — 2/Ak' as (52,3/|S^2| ^ 1- There- 
fore, it is not possible in conventional nonlinear optics to 
get the complete conversion, 773 ^ 7710, without careful 
phase matching. 

If such a compensation is performed, i.e. if 
(K/2Afc') > 1 (by proper choice of detunings 52 and (S3), 
optimum conversion occurs for a density-length product 



I IT (5?<S3 
Nz\ ^ -m = TT ^ , 

i°pt 2 7MT7^|f]2| 



At this optimum propagation distance, the relative ab- 
sorption losses are given by the parameter F, Eq. (pO|). 



Thus, we conclude that in the limit of undepleted drive 
field(s) both EIT and maximum coherence schemes per- 
form much better than conventional nonlinear optics, re- 
garding the number of atoms -^-s^j^p^ necessary for opti- 
mum conversion and robustness to phase mismatch. The 
aim of the present paper is to investigate whether the em 
energy can be transferred to the generated wave from 
both pump fields, and whether the attractive features of 
EIT- and maximum coherence-assisted conversion sur- 
vive when the drive field is depleted. This requires the 
solution of the complete nonlinear propagation problem 
which we treat with the Hamiltonian formalism outlined 
in the following section. 



III. HAMILTONIAN APPROACH 

The resonant interaction of em fields with atomic sys- 
tems is described by Maxwell's equations for the fields 
and a set of master equations for the density matrix of 
the atoms. The solution of this coupled set is rather dif- 
ficult and except for some very special cases impossible 
analytically. In order to derive analytical solutions for the 
field propagation, approximations are needed that allow 
to eliminate the atomic degrees of freedom and to express 
the polarization in terms of field variables. A common 
approximation, which is not perturbative in the fields, is 
the adiabatic solution of the density matrix equations. 
This straightforward approach provides, when applica- 
ble, full information about both atoms and fields in the 
adiabatic limit. 



A. General formalism 

Instead of the usual explicit scheme of adiabatic elim- 
ination which is rather cumbersome, we use here a dif- 
ferent, implicit approach ^ which yields directly 
effective field equations. This approach is based on the 
representation of the medium polarization P as a partial 
derivative of the time-averaged free energy density of a 
dielectric with respect to the electric field E : 



P = -( N 



OH 
'dE 



(33) 



where (...) denotes quantum-mechanical averaging, and 
H is the single-atom interaction Hamiltonian. For the 
field given by Eq. (^, we can write: 

(dH \ 
^^d£^ exp(-z(^,t - k,z)) + c.c.j , 



(32) so that the propagation equation becomes: 



d£i „ (^i / dH 



dz 



d£* 



(34) 
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We here consider light-atom interaction processes that 
are adiabatic, that is the atomic system can be assumed 
to follow the evolution of the instantaneous eigenstates. 
For example, if the atomic system is at some initial time 
^0 in the nondcgcncrate eigenstate \ipo{ta)) of the inter- 
action Hamiltonian, i.e. 



(35) 



which is at t = to usually identical to the ground state of 
the atoms, it will remain in this state at all times. Fur- 
thermore, we disregard irreversible dissipation processes. 
As we show later, for the processes we consider here, this 
is justified even for rather long pulses, i.e., longer than 
the natural lifetime of the excited state |3). In this 
limit one finds 

/dH\ dH dXo 

Hence the propagation equation can be written as: 



9£ 

dz 



■' = —i2tt — -N- 



d£* 



(36) 



For the following it is useful to express the field amplitude 
£j in terms of photon flux ijj , Eq. ( [l^ ) , and phase ipj : 

£j = \£j\exp{-ipj}. 

Separating the real and imaginary parts, we find from 
Eq. (H): 



(37) 



dr]j _ m' 
dz difj ' 

difij _ dW 
dz drjj 



These equations have the form of Hamilton equations of 
classical canonical mechanics with action and angle vari- 
ables rjj, and ipj, "time" z, and the Hamiltonian function 



B. Constants of motion for resonant 3- wave mixing 

For atomic systems with a closed loop of transitions, 
the set of canonical equations ( p7[ ) can be further sim- 
plified and in fact under some conditions explicitely in- 
tegrated. This is the case when the number of coherent 
fields involved is not too large. To illustrate the proce- 
dure, let us consider the eigenvalue equation (|3^) for the 
three- level system in Fig. 1. In rotating- wave approxi- 
mation, the light-atom interaction Hamiltonian is given 
by: 



H 



S2\2) {2\+S,\3) (3| 



(38) 



and the eigenvalues are determined by the characteristic 
equation: 

Ao (^2 + Ao) (^3 + Ao) - {nl + + nl) Ao - 

-nl53-filS2 = -2nin2i^3coaip. (39) 

Here fi^ are the Rabi frequencies, related to the photon 
flux via the coefficients ^j, Eq. (p^), as Qj = ^/JlfffJ. 
The relative phase ip of the em waves is 



(40) 



which includes the residual phase mismatch Ak. Also 
the multiphoton resonance condition 



UJ3 = UJl + UJ2 



(41) 



has been used. 

One can see from Eq. ( p9| ) that Aq and, hence 7i', 
depend on the field phases ipj only through the relative 
phase ip. Therefore, we have: 

difi dif2 d(p3 V dip J ' 

An immediate consequence of this symmetry of Ti.' is the 
existence of constants of motion. Substituting the above 
equations ( p2| ) into the first line of Eqs. (^) yields the 
well-known Manley-Rowe relations M : 



drji _ dr]2 _ 9??3 
dz dz dz ' 



(43) 



which correspond to two independent constants of mo- 
tion: 



V1-V2 = Vio - V20- 



(44) 



Here rjjQ = rjj^z = 0) are the photon flux values at the 
entrance to the medium. Taking into account the multi- 
photon resonance condition (p4|), one finds furthermore 
that the total intensity of the em fields is conserved: 
Ii + 12 + 13 = const{z). The Manley-Rowe relations and 
the constants of motion tell us that in the process un- 
der consideration, the energy is transferred from the fre- 
quency components wi, 0^2 into and back, with equal 
rates and without losses, which of course is expected for 
a dissipationless nonlinear medium. 

The relations Eq. ( ^ ) enable us to re-write r]j as: 



V2{z) = V2Q - J{z), 

m{z) = 7730 + J{z). 



(45) 



-nni |1) (2| + ni72e"^ |2) (3| + nQ.3 |1) (3| -I- H.C., 



J{z) characterizes the amount of energy exchange be- 
tween the waves and has the initial condition J{z = 0) = 
0. 
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Thus the original problem with six amplitude and 
phase variables can be reduced to two variables J and 
(ys by a canonical transformation. This leads to 



dJ_ 

dz 
dip 
dz 



m 

dip ' 

m 
a/' 



with the new Hamiltonian function 



AfcJ= -NX. 
2 



(46) 
(47) 

(48) 



As can be seen from Eqs. ( p^ ) and (|39|), Ti. (or A) does 
not depend on the coordinate z explicitly. Therefore, 
TL (or A) is a fourth constant of motion expressing the 
conservation of the energy density of the medium with 
respect to z. 



C. Solution of the wave- propagation problem 

To solve the remaining two equations of motion for 
J{z) and ip{z)^ the Rabi-frequencies ilj are expressed in 
terms of rjjQ and J, and the characteristic equation (|39| ) 
is written in the form 



G(A, J) — g[J) cos (/5. 

Differentiating both sides with respect to tp yields 

dG dGdX . , — — 

d^ = ^d^ = -''''''^ = ^^'--^ ■ 

Substituting this relation into Eq. (|4q), we find: 



(49) 



dz ~ 2 dG/dX ' 



(50) 



The choice of the sign in Eq. ( pO| ) depends on the sign of 
sin tf a,t z — 0. Integration of Eq. (^0|) gives an implicit 
solution for J{z): 



_^7V^_ j dG{J') dJ' 



(51) 



To analytically evaluate the remaining integral, we note 
that both functions — G^ and dG/dX are polynomials 
in J: 



g = -2^/Xi/Li2At3A/(?/10 - J) (??20 - J) (??30 + J), 



G ^ Go 
1 



— - V a J" 

dx " ^ ™ 



771 — 



Therefore, equation (|5C|) describes a one-dimensional fi- 
nite motion of a pendulum in an external potential. The 



solution is in general given by some combination of ellip- 
tic functions |18| with parameters determined mainly by 
the roots J„ of the polynomial equation: 



(J) - (J) = 0. 



(52) 



The allowed range of J, corresponding to the region of 
classically allowed motion of the pendulum, lies between 
zero and the smallest positive root J\ of the polynomial 



The eigenvalue A is a constant of motion (cf. Eq. ([48D), 
and can thus be found from the characteristic equation 
( ^ ) with parameters taken at the medium entrance z — 
0: 



Go (A) = g{z = 0) cos <p{z = 0). 



(53) 



Thus, we have reduced the propagation problem to 
solving two algebraic equations: ( p3| ) for A and (^2|) for 
the roots J^. If this can be done explicitly, the Hamilto- 
nian method provides an analytical solution to the prop- 
agation problem. But even if an explicit solution is not 
possible, it considerably simplifies numerical calculations. 
Apart from the advantage of being convenient from the 
formal point of view, the Hamiltonian method allows a 
deeper insight into the nature of the wave propagation 
process. For example, the method provides a direct ac- 
cess to the stability analysis of the solution by referring 
to a well-developed theory of Hamiltonian systems . 

To understand the dynamics of the system, in particu- 
lar interesting quantities such as the conversion efficiency, 
it is useful to discuss the physical meaning of the coefh- 
cients Am and Om- Considering the canonical equation 
( ^ ) for the relative phase: 



dp _N dX _N dG/dJ _ N Ai + 2^2 J + 3^3 
97 ^ Y9J ^ Y dG/dX ^ y ao + aiJ + a2J2 



(54) 



one recognizes that the Am and am describe the linear 
and nonlinear refraction coefficients of the medium. E.g. 
if J is sufficiently small, the first term NAi/2ao on the 
right-hand side of Eq. (|5j) can be identified with the 
phase mismatch induced by the linear refraction, includ- 
ing both contributions from the three-level interaction 
and the residual mismatch Afc. 



D. Generation of field with frequency ujs 

In the context of the present discussion we are most 
interested in the generation of the uj^ mode from vac- 
uum 7730 = 0. In this case, the eigenvalue equation (|5^ ) 
reduces to 



Go(A) = X{X + 62)iX + S3) - ^ilr|loiX + S3) 

-^2?/2oA = 0. 



(55) 
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The nonvanishing coefRcients in the expansion of G and 
dG/dX are given by 

Ai = q\{X + 52){X + 53) + XiX + S2) + XiX + S3) 



-Mi'7io - M2?/2o (56) 
(A + S3) + 1.12X - (A + ^2) , 
{3X + S2 + S3) + q {fii + fi2 - f^s) , (57) 



A2 
A3 
ao 
ai 
a2 

where 



= q', (58) 

= 3a2 + 2A {S2 + S3) + S2S3 - ^llmo - ^^2V2o, (59) 

= 2q{'sx + S2 + S3) + i^il + ^i2-^l3), (60) 

= 3g^ (61) 



q = -2Ak/N. 



(62) 



With this one finds for the denominator in the integral 
(pT|), which determines the allowed values of J: 

-G"^ = 4:fJ,ifi2f^3J ima - J) {V2a - J) 

^{Ai+A2J + A3J^)^J^^0. (63) 

The second term in this expression is non-positive, so the 
smallest positive root of the polynomial is bounded by the 
minimum of 7710 and 7720- This reflects the fact that the 
conversion process stops when the energy of the weaker 
of the two pump fields is entirely depleted. In order to 
reach this limit and thus in order to attain maximum 
conversion efflciency, the second term in ( |6^ should be 
small, which corresponds to the phase mismatch to be 
negligible: 



A1+A2J 







(64) 



In order to see what parameters are required to approx- 
imately satisfy this condition, we have to analyze the 
coefficients Am- 

It is worthwhile to point out that the behavior of the 
total phase is described by Eq. (Eol): 



(65) 



, , G AiJ + A2J^ + A3J^ 
cos (p(z) = — = — . 

9 2^y^llfl2^J■3J ivw - J) imo - J) 



When the total phase mismatch is compensated, Ai + 
A2J + A3 ~ 0, the phase follows the equation: 

cos ip{z) = 0, 

that is, the phase is constant and equal to +n/2 or 
— 7r/2 when all of the waves are present, and jumps by 
TT to — 7r/2 or -t-7r/2 when one of the wave intensities ap- 
proaches zero. 

As it is seen from Eqs. ([57|), (^), the terms A2, A3 
responsible for the intensity-dependent refractive index 
are both proportional to the residual phase mismatch 
Afc. One can show that in atomic (molecular) media, 
where the background refractive index is given by Eq. 



(p5|), Afc is always sufficiently small so that the lead- 
ing term in the phase mismatch is Ai, while A2J and 
^3 are negligibly small. It is also true, that the influ- 
ence of the intensity-dependent phase mismatch remains 
insignificant even when the linear refraction is compen- 
sated: Ai = 0. Therefore, in what follows, we will disre- 
gard the terms A2J and A3 J^. In this case, the solution 
of the propagation equation (^ gives the following de- 
pendence of J{z) in implicit form: 



±KZ + xo = F['yi{J),p] 



ao 



{F [12 {J), p] 



+ (l-f)^[72(J),P^p]}, 



(66) 



where xo is an integration constant, and F{'^,p) and 
n(7,r,p) are the elliptic integrals of the first and third 
kind, respectively [ |l8| . k is the nonlinear conversion co- 
efficient defined as 



N ^/^ll^l2^J.zJ2 

K = , 

2 ao 

The parameters of the elliptic integrals are: 



71 (J) = arcsin-i/ — , 



(67) 



lj2{Jl-J) 
Jl {J2 - J) ' 



72 (^) = arcsin . 



with Jl and J2 being the roots of Eq. (r3 



(68) 



J2^ 



Bi 



1 , 

2 (^10 + '720 + Bi) 



^^yivw + mo+Bi) -477107720, (69) 
AI 



4^1/^2M3 

In the next section we discuss the solution of the prop- 
agation problem for several field-atom interaction config- 
urations resulting in the EIT-assisted, maximum coher- 
ence, and conventional nonlinear-optical regimes of the 
frequency conversion. 



IV. ANALYTICAL SOLUTIONS FOR 
RESONANT 3- WAVE MIXING 

We assume in our simple 3-level model system that 
the transitions |1) — |3), |2) — |3) as well as |1) — |2) are 
allowed. To take into account the much weaker coupling 
between levels |1) and |2) which in reality is a two-photon 
transition, we assume <^ 112- An analytic solution 
for the more general 4-wave mixing system can also be 



If 



The 



obtained within the Hamiltonian formalism 
explicit results are however rather lengthy and not very 
instructive. 
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A. EIT-based upconversion 

The EIT scheme imphes two-photon resonance be- 
tween the states |1) and |3): = 0. According to the 
idea of EIT-assisted nonlinear optics, we suppose here 
that r]2o > rjiQ. Since /X2 ^ /Xi, we always have ^ ^^i- 
We allow, for the moment, for spontaneous decay from 
state 1 3) out of the system, but require the initial Rabi 
frequency $720 to be much larger than the decay rate: 
^20 ^ 7- Under these conditions, the eigenstate that 
asymptotically connects to the ground state of atoms for 
t ^ — oo is the state 



lV'o(^ = 0)) 



|i> 



■lip 



10 

corresponding to the complex eigenvalue 



"20 



|3> 



A: 



"lO 



Since flw ^ one has: |A| <C 7. Correspondingly, the 
probability of spontaneous decay of the state lipo), which 
can be defined as Im(A)T with t being the characteristic 
pulse length, remains small even for rather long pulses: 
T ^ To = (QIq/QIq) "f^^ with To > 7"^. For such 
pulses, we can safely disregard the spontaneous emission 
and put A « 0, which means that the corresponding adia- 
batic state is the usual dark state of EIT. Moreover, since 
flio ^ ^20 the eigenstate lil^oiz = 0)) is an only slightly 
disturbed ground state |1). Therefore, the requirement 
of adiabaticity of the process is satisfied automatically at 
the medium entrance, z = 0. 

The coefficients Ai and a™ are given for the regime of 
EIT by the following expressions: 



ao — -M2??20, 

ai = 2qS2 + (/ii + /i2 



Ms) ■ 



(70) 
(71) 
(72) 



1. Undepleted EIT-generating field 



In the original proposal of EIT-assisted frequency con- 
version, the UJ2 field is assumed to be very strong and 
undepleted. The latter condition implies 7710 <C 7720- In 
this situation the solution ( |66| ) can be well approximated 
by 



mo 



1 



Afc' 
2k, 




N 

2 V /i2'?20 



Ak' 



Ak 



N 11362 

2 M2?720 ' 



,(73) 

(74) 
(75) 



which, of course, coincides with that obtained in Sect. 
II, Eqs. (|lj) and (|l5|), under the undepleted drive ap- 
proximation. The influence of the intensity-dependent 
phase mismatch (terms A2J and A^J^) is negligible in 
this case, as was discussed above. 

One can immediately see from Eqs. ([70|), (|57|), ( |5^ ) 
that for vanishing residual phase mismatch q = (which 
can be realized by adding a buffer gas with proper disper- 
sion and for 62 — there is perfect phase matching: 
Ai = A2 = A3 = 0. Consequently maximum energy 
transfer into the generated wave is possible. This 
applies also for the case when linear refraction is com- 
pensated by a detuning 62'- 



52 = 



2 /i2'720 



Afc. 



(76) 



N fl3 

For both of these phase matching procedures, the max- 
imum value of 773 = J achieved at Ze = ■n/2Ke is equal 
to ryio. Thus the instantaneous fractional conversion ef- 
ficiency, defined as 



min (7710,7720) 



(77) 



can become unity. It should be noted, that e = 1 is in 
general only achieved at a specific instant of time since 
Kg depends on time via 7720 {t — z/ c) = 7720 /20 {t ~ z/c) 
(7720 is the amplitude and /20 {t — z/c) is the temporal 
envelope of the 7720 {t — z/c) pulse). Hence the transfer 
of energy is complete only for the part of the wi pulse 
which satisfies the condition ■\/ /20 ~ z/c) « 1. A high 
overall conversion requires a flat temporal profile of the 
"coupling" pulse aj2- This is easy to implement, however, 
since always 7/3 = J <C 772 (~ 7/20) and 771 <C 7720- We have 
then f2i, r23 <C Q.2 for the whole medium which means, in 
particular, that all the atoms throughout the medium are 
essentially in the ground state. Therefore, the 0^2 pulse 
can be arbitrarily long. 

In contrast to the fractional conversion efficiency, Eq. 
([77]), the total conversion efficiency, defined as 



W{z) = 



At L03rj3{z,t) 



dt ujirii{z,t) + ijj2ri2{z,t) 



(78) 



remains very small because 7710 <C 7^20- 



2. Depletion of the pump fields 

Conversion with e sa 1 can be achieved, however, also 
for the case when both pump fields uji and W2 are of com- 
parable intensity. For example, for exact phase matching 
and 7710 = 7^20 = 7^0 (which would correspond, e.g., to the 
second harmonic generation), the solution Eq. ( |66[ ) can 
be written as: 



M2 

— / 
Ms 



Arth 




(Mi + M2 - M3 ) 



J 



(79) 
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which demonstrates that J monotonically approaches rjo 
as z increases. The form of this solution for given re- 
tarded time is shown in Fig. |[ The monotonic de- 
pendence ensures that after sufficiently long propagation 
length all parts of the pulse 7710 — 7720 Sire converted into 
the 073 wave so that not only e « 1, but also the total 
energy conversion efficiency W reaches unity. 




FIG. 2; EIT-upconversion for perfect phase matching. Solid 
line: spatial evolution of J according to Eq. (^9|) for given 
retarded time. For sufficiently large propagation length there 
is complete conversion. Dotted lines: corresponding spatial 
evolution of the atomic bare state probabilities. The ampli- 
tude of the excited state |3) is negligibly small (not shown). 
Parameters: /i2/M3 = 0.5, /ii//^3 = 0.05. 

When the phase mismatch is not compensated, the 
solution for 7710 = 7720 = Vo is given by more gen- 
eral expression Eq. ( |66|) with Bi = rjo (AA:'/2Ke)^ and 
J2.1 ~ % (1 i Afc'/2Ke). The spatial dependence of the 
generated field intensity (Fig. ^ differs only slightly from 



J - 7/0 1 



Afc' 



1 



Ak' 



where sn[x;p] is the Jacobi elliptic sine function. The 
period of intensity oscillations seen in Fig. ^ can be esti- 
mated as 18 



very 



Ze ~ ln(16Ke/Afc') . 

Since in EIT regime we have always: Afc' <C 2Ke 
high total conversion efficiency is achieved without phase 
matching also for 7710 ~ 7720 when the drive fields are 
substantially depleted. 

3. Evolution of the atomic state 

In the case of EIT with comparable intensities of the wi 
and UJ2 pump fields, the intensity of the wave increases 




FIG. 3: EIT-upconversion with finite phase mismatch. Solid 
line: spatial evolution of J. Dotted lines: corresponding spa- 
tial evolution of the atomic bare state probabilities. The am- 
plitude of the excited state |3) is negligibly small (not shown). 
Parameters are the same as in Fig. Q, Afc'/2Ke — 0.06. 



with z and eventually exceeds the one of the wi, UJ2 waves. 
Therefore, the adiabatic state of the atomic system will 
no longer approximately coincide with the ground state 
|1) over the entire medium length at a given retarded 
time. It will evolve via different superpositions of the 
bare states |1) and |2) following the change of the fields 
(Fig. I and |) [|o|- We see that the "conventional" 
EIT regime (one strong and one weak field, and thus one 
highly populated and one slightly populated state) does 
not hold over the full conversion cycle. We should stress, 
however, that the atoms remain in the adiabatic state 
corresponding to A = 0, i.e., in the dark state. There- 
fore, all the features of the EIT-assisted conversion pro- 
cess such as vanishing linear absorption and refraction 
as well as resonantly enhanced nonlincarity are present 
even when the pump field is depleted. 

When most of the input energy is transferred to the 
^3 pulse, i.e. for larger propagation distances into the 
medium, the dark state It/jo) approaches the bare state 
|2) 

iV'o) |2), 

which does not coincide with the initial atomic state (i.e. 
the ground state |1)). Thus we encounter two problems: 
First, since all atoms are in the ground state before the 
arrival of the pulses, there must be an adiabatic transfer 
in the initial phase from |1) to IV'o)- This requires a spe- 
cific time-order of the pulses. In particular VI2 should ar- 
rive after the other pulses. This time ordering is however 
not reflected in our solutions. Second, after the inter- 
action with the pulses, there is energy left in the atomic 
system. This is in apparent contradiction to the fact that 
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full conversion of photon energy has taken place and no 
energy got lost. It should be noted though that the size of 
the effect is rather small in the considered limit of a large 
initial number of photons. This contradictory behavior is 
a result of the limited validity of the assumptions made. 
We have here completely neglected non-adiabatic correc- 
tions, which are small on the level of individual atoms, 
but add up in the field evolution when integrated over 
the whole sample. The most important effect ignored by 
this is the group velocity delay of the weak fields. 

Although the group delay cannot be treated within the 
current approach, it may resolve the problems. In the 
first part of the medium, i.e. for small Nz the generated 
field is still small. Thus the group velocities of the pulses 
fii and r^s should be smaller than that of the strong 
drive field VI2 ■ Then the leading edge of the ri2-pulse will 
always propagate ahead of the other pulses, guarantee- 
ing that the adiabatic state |V'o) asymptotically matches 
with the ground state at early times for all values of z. 
For larger propagation distances into the medium and 
for times in the central part of the pulses a substantial 
portion of the energy is transferred to the field fia. Now 
the situation is reversed. fJs is strong and f22 is weak. 
Thus the pulse ^2 will have a smaller group velocity and 
should lag behind at the trailing edge of the pulses. Con- 
sequently also for large times the adiabatic eigenstate 
|'0o) approaches asymptotically the ground state |1). 

A full quantitative account of the leading order non- 
adiabatic corrections requires however a reformulation of 
the Hamiltonian approach. This is beyond the scope of 
the present work and will be presented in a future publi- 
cation. 



B. Maximum coherence case 



medium, z = 0, are: 



Cl 



C2 = 



10 



A 



VA^TTI?;' 

C3 < Ci,C2. 



(84) 
(85) 



Similar to the EIT case, the probability of spontaneous 
decay of the state |-0o) is small: Im(A)T <C 1 for large 
detuning assumed here, 83, 3> 7, even for pulses that are 
longer than 7^^: t <C ('5|/ri2o) 7^^- Again, we will con- 
sider only such pulses and neglect spontaneous emission. 

The situation of the maximum coherence, |ci| = |c2| = 
l/\/2, takes place when the Rabi frequency exceeds the 
detuning on transition |1) — |2) including the ac Stark 
shift induced by the a'2 field: fiip ^ |(52 — il2o/'^3|- 
In this case A « ilio and the atomic state at the 
medium entrance is the superposition: |V'o(^ = 0)) = 
(|1) — 1 2)) 1^/2. However, in order for the atoms to get 
there from the ground state, one first has to have the 
inverse situation 82 ^ i^io, and then adiabatically de- 
crease 82, simultaneously increasing This can be 
accomplished by, e.g.. Stark-chirped rapid adiabatic pas- 
sage |8|. We will not include this process explicitly in 
the consideration, but assume that when the L02 pulse 
arrives, the condition fiio S> 82 is already satisfied. As 
a consequence, the energy taken from the leading edge 
of the LOi pulse for the superposition preparation, is not 
taken into account. However, if the number of photons 
in the oji pulse is much larger than the total number of 
atoms along the field propagation path, these preparation 
energy losses (?ia;i/2 per atom) can be neglected. 



The mechanism of the nonlinear optics with maximum 
coherence is realized when the detuning 83, is much larger 
than the other parameters including all Rabi frequencies 
and 82- In this case the coefficients Ai and are as 
follows: 



Ai ~ q[{2\ + 82) 8s. ~ ^Jilr|w - ^l2r|2^^\ 

+fJ'i83 + M2A - ^3 (A + 82) , 
qq ~ (2A + (^2) (53 - Mi?no - Ai2?720, 
ai = 2q83 + (/ii + Ai2 - M3) • 

1. Preparation of maximum atomic coherence 
The energy eigenvalue for large 83 is given by: 



(80) 

(81) 
(82) 




"20 



where 8^ = 83 + 1"/. The amplitudes of the adiabatic state 
corresponding to this eigenvalue at the entrance to the 



2. Undepleted coherence- generattng field 

When a strong, undepleted uoi field is supposed: 7720 ^ 
Tyio (implying also undepleted coherence on |1) — |2) tran- 
sition), the generated field intensity experiences sinu- 
soidal oscillations with respect to the propagation length: 



J = 



mo 



■ sm 



l^m. ^ A / 1 



2Km 



A/c' 




Ak' = 



in correspondence with the analysis of Sec. 

|C2|'= 



(86) 

(87) 
(88) 



II, Eqs. (y), 

P12 = 1/2. 



( P2| ) and (|24|), and assuming |ci| 

Since Km is of the order of Afc', substantial transfer 
of energy from the 1^2 field into the W3 field occurs even 
without compensation of the phase mismatch. If one 
can manage to make Ak' small: Ak' <ti k,„ then the 
best possible regime for conversion can be attained: not 
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only complete transfer for the pulse maximum is realized: 
Ji « Jmax = '720 j but also the conversion length I — 
does not depend on time, and hence, all parts of the uj2 
pulse are homogeneously converted into the UJ3 wave. 

There are several possibilities to reach this regime. 
First, one may set Ak to zero by use of the buffer gas, 
and choose the input laser parameters such that: 



(M2 - Ma) 
Ml 



(89) 



This condition can be satisfied only for the time inter- 
val when •\//io {t — z/c) « 1. That is, a flat temporal 
profile of the oji pulse is required. Since this is not al- 
ways possible, another method can be used. One may 
simply increase the intensity 7710 so that the refraction 
contribution of the |1) — |2) transition becomes small: 
\/ Mi/'7io *C (M2 ^ Ms) /'^3i ^-i^d choose the detuning ^3 as 



(M2 - Ma) 
2q 



(90) 



Since in this regime rj^ = J 7710, we have at any 
propagation distance: fli » ^2/^3, ^l/Ss- There- 
fore, the atoms remain in a superposition \ijjQ) — 
(|1) — 1 2)) /\/2 throughout the medium. 

It should be noted that although the instantaneous 
conversion efficiency e can reach unity, the total energy 
conversion efficiency is always small, W 1, because 
Ji « V20 < mo- 



3. Depletion of the pump fields 

The situation is completely different for comparable 
intensities of the two pump waves: 7720 ~ ??io • 

For the maximum coherence regime, |ci| = |c2| = 
l/-\/2, to be attained, the condition Q.i ^ must 
be satisfied (see above), which reduces to ^/i2?/io ^ 
<^3\/mi/M2 at 7720 ~ ^10- That is, the intensity of the 
pump pulses cannot be large, and the refraction contri- 
bution of the |1) — 12) transition (^ Mi/^10) is not small, 
so that AA:' « [N/4) \J pi\lr\\Q >• Km. Thus, the "exact 
maximum coherence" regime fails when the drive field is 
depleted. 

However, the phase matching can still be performed 
for 7720 ~ r\\Q if the requirement for |ci| = |c2| = l/\/2 
is lifted. One of the possibilities is to cancel the resid- 
ual refraction by addition of buffer gas: A/c = (then 
^2 = ^3 = 0), and to choose the interaction parameters 
such that A\ w 0. The latter condition can be satis- 
fied for specific relations between the detuning 82 and/or 
(^3, and the input intensities 7710, ?y20- For example, for 
^10 = ^20 = % a-iid when /i2 < ^^^^ one may choose the 
detunings ^2 = and 83 such that 



M3??o = 



Ml 



(M3 - M2) 



(91) 



Taking into account condition (^l|), the eigenenergy A is 
reduced in this regime to: 



A = VMi% 



M3 



(M3 - M2) 



(92) 



and the corresponding adiabatic state at the entrance to 
the medium is: 



- 0)) 



1 



|1>-C|2)). (93) 



Thus, the atoms are prepared not with maximum co- 
herence, P12 = 1/2, but in a specific superposition with 
P12 = C/(l + C^) which is also not small. 

The propagation problem solution in this case resem- 
bles that of the EIT regime, Eq. (^): 



2 



c 



Arth 




J \ _ Ml + M2 - M 3 / J 

M3 



M3 



(M3 - M2)' 



(94) 



(95) 



which is demonstrated in Fig. ^. The conversion 
coefficient differs from that of the undcpleted coher- 
ence case Km, Eq. (|87|), only by a numerical factor 
2v/(M3 -M2)/M3- 




FIG. 4: Upconversion with initial maximum coherence. Solid 
line: spatial evolution of J according to Eq. (^^ for given 
retarded time. For sufficiently large propagation length there 
is complete conversion. Dotted lines: corresponding spatial 
evolution of the atomic bare state probabilities. The ampli- 
tude of the excited state |3) is negligibly small (not shown). 
Parameters: /i2//i3 = 0.5, /ii//i3 = 0.05. 

As we can see, although the total conversion efficiency 
may approach unity, the robustness of the maximum co- 
herence scheme to the phase mismatch is lost. One has 
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to fulfill the phase matching conditions (like (|9T|)) with 
sufficient accuracy to get a large conversion. A consid- 
erable transfer of energy into the generated wave can be 
expected when the envelope of the pump pulse satisfies 
the condition: 



1 - Vfio {t - z/c) < 



(M3 -Ai2)' 



and the detuning S3 
given by condition (te 



SI 



5'o differs from the value So, 



by not more than 



(M3 



AO 



Both of these conditions require that the values of the 
transition coupling constants, fi2 and fi^, be close. 

Similar to the EIT regime, the adiabatic superposi- 
tion does not coincide with the state \ipo{z = 0)), Eq. 
(p3|), over the entire medium, but will evolve in the 
course of propagation, following the change of the fields 
(Fig. 1^). For large density-length products Nz when 
fii <C ^3/ S3, the atomic state tends to the ground state: 
IV'o) ~* Thus, "the maximum coherence" is not 

maintained throughout the light propagation path, al- 
though it is still sufficiently large for most part of the 
medium. One can see from Fig. H that as the coherence 
gets smaller, the conversion slows down, and for very 
large Nz (hence, for small atomic coherence), the gener- 
ated intensity only slowly approaches its maximum. 

Here again we encounter an apparent contradiction re- 
lated to the limited validity of the adiabatic approxima- 
tion. The atoms do not return to the initial state despite 
the overall energy conservation. Now the preparation en- 
ergy for the atomic coherence is extracted as well. Under 
the conditions discussed here, i.e. large photon number 
as compared to number of atoms, this is however a small 
effect. A resolution of this problem should be found when 
non-adiabatic corrections are taken into account, which 
will be the subject of further investigations. 



C. Conventional nonlinear optics (weak excitation) 

Finally, we consider the regime corresponding to con- 
ventional nonlinear optics. This regime takes place at 
weak excitation when both detunings S2 and ^3 are 
much larger than all Rabi frequencies. Then we have 
A iliQ/S2 and \iPq) w |1) for all atoms in the medium. 

The coefficients Ai and Um do not depend on laser 
intensities in this case: 



Ai ~ qS2S3 + H1S3 - ^3^2, (96) 

ao ~ 5263, (97) 

ai = 2q(52 + '53) + (/^i-f M2-M3)- (98) 

Since ai min (?7io, ?72o) ^ cto, the solution Eq. (pq) can 



be well approximated by 



J = Jisn^ 



N V/^l/i2M3 J2 

2 S2S3 



(99) 
(100) 



with J2,i determined from Eqs. {p9\j, (|96|). This formula 
is valid for any relation between 7710 and 7720. 

When the phase mismatch is compensated: Ai « 0, 
which is made by proper tuning: 



("3 Ml 

03 02 



(101) 



the roots J2,i are: Ji = min (7710, 7720), J2 = 
max (7710, 7^20) ■ Therefore, the optimum conversion can 
be realized also in the case of the off-resonant nonlinear 
optics. In particular, for 7710 « 7720 the complete transfer 
of energy to the generated wave occurs: 

J = 77osn^ [Knz;p ^ 1] « 770 tanh^ (k„z) , 

leading to perfect total conversion efficiency, W ^ 1. 
However, the nonlinear conversion coefficient Kn is much 
smaller in this case than in the EIT and maximum co- 
herence regimes (Kg and Km, Eqs. (^ and (|87|), respec- 
tively). Therefore, the energy transfer is accomplished 
for much larger density- length products Nz. 

At last, we note that if the phase mismatch is not com- 
pensated, we have Bi ^ (t^iq -I- 7720) and the solution is 
reduced to traditional formula of nonlinear optics |B|: 



J = N 



MiM2M3'7io??20 sin^ (Afc'z/2) 



AS^Si 



(Afc'/2)' 



(102) 



with the total phase mismatch Ak' = Ak — N ^1/282 
Nfi3/253. 



V. SUMMARY 

Resonant optical processes based on atomic coherence 
effects such as EIT or maximum coherence allow for a 
maximum nonlinear conversion of photons within a much 
smaller density-length product than possible in schemes 
of conventional off- resonant nonlinear optics. Since de- 
grading mechanisms such as phase-mismatch and ab- 
sorption scale with the same density length product, 
EIT or maximum coherence can thus substantially re- 
duce the requirements for efficient nonlinear optics. As 
a measure for the efficiency of using atoms for nonlin- 
ear conversion processes, one can introduce the ratio 
of the number of converted photons n ^ ATjfjjo, to 
the number of atoms Nat = NAl ^ NA/k needed {A 
is a cross-sectional area and Tj is the duration of the 
pulse). This figure of merit n/Nat is largest for the EIT 
scheme: n/Nat = ''"1^10, considerably smaller for the 
maximum coherence method: n/Nat — T2Q20 (^2o/<53) 
(with ^20/63 1), and is very small for conventional 
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nonlinear optics: n/Nat = ^2^20 {^10^20 / 5^52) (with 
^2Qlh < 1 and f^io/<52 < !)■ 

In previous theoretical studies of EIT- and maximum- 
coherence based nonlinear optics, undepleted drive fields 
or a constant coherence where assumed. If one takes into 
account the resources to maintain the drive field or the 
constant coherence, the overall efficiency of the processes 
is in this limit tiny. We have shown in the present pa- 
per that when this restriction is lifted and coherence- or 
transparency maintaining fields with comparable inten- 
sities to the pump fields are considered, it is possible to 
achieve also a maximum overall conversion. To study the 
conditions for this, we have derived analytical solutions 
for the pulse interaction in the adiabatic limit using an 
approach that maps the propagation problem to that of a 
nonlinear pendulum. Under certain conditions this prob- 
lem could be explicitly integrated allowing for a simple 
discussion of the physical processes involved. 

We have found that when the coherence- or trans- 
parency maintaining fields are depleted, the atomic state 
does not remain constant but evolves along the propaga- 
tion path following the change of the fields. Therefore, 
the conversion process cannot be associated with maxi- 



mum coherence or EIT in traditional sense alone. Never- 
theless, the main features of EIT, - reduced linear absorp- 
tion and refraction, and enhanced nonlinearity are still 
present in the latter mechanism. Although in the " maxi- 
mum coherence" mechanism the robustness to the phase 
mismatch is lost, the nonlinear conversion coefhcient re- 
mains much larger than in conventional nonlinear optics. 
Thus, for both mechanisms, a complete conversion can 
be achieved within a small density-length product. 

We have also encountered several limitations of the 
strong adiabatic assumption used in the Hamiltonian ap- 
proach. An extension of the approach to take into ac- 
count non-adiabatic corrections and thus such important 
effects as group delay is currently under investigation and 
corresponding results will be presented elsewhere. 

The work of E.A. K. was supported by Alexander 
von Humboldt Foundation. We would like to thank R. 
Unanyan for bringing to our attention the Hamiltonian 
formalism in nonlinear optics and the relevant publica- 
tions jl^, |l2|, and him and K. Bergmann for many use- 
ful discussions. E.A. K. would also like to thank K. 
Bergmann for the hospitality and his continuous support. 



[1] S.E. Harris, Physics Today 50, No. 7, 36 (1997). 

[2] S.E. Harris, J.E. Field, and A. Imamoglu, Phys. Rev. 
Lett. 64, 1107 (1990). 

[3] M. Jain, H. Xia, G.Y. Yin, A.J. Merriam, and S.E. Har- 
ris, Phys. Rev. Lett. 77, 4326 (1996); S.E. Harris, G.Y. 
Yin, M. Jain, H. Xia, and A.J. Merriam, Phil. Trans. R. 
Soc. Lend. A 355, 2291 (1997). 

[4] M. Jain, G.Y. Yin, J.E. Field, and S.E. Harris, Opt. Lett. 
18, 998 (1993). 

[5] G.Z. Zhang, K. Hakuta, and B.P. Stoicheff, Phys. Rev. 
Lett. 71, 3099 (1993); G.Z. Zhang, M. Katsuragawa, K. 
Hakuta, R.I. Thompson, and B.P. Stoicheff, Phys. Rev. 
A 52, 1584 (1995); G.Z. Zhang, D.W. Tokaryk, B.P. Sto- 
icheff, and K. Hakuta, Phys. Rev. A 56, 813-819 (1997). 

[6] C. Dorman and J. P. Marangos, Phys. Rev. A 58, 4121 
(1998); C. Dorman, I. Kucukkara, and J. P. Marangos, 
Phys. Rev. A 61, 013802 (1999). 

[7] A.J. Merriam, S.J. Sharpe, H. Xia, D. Manuszak, G.Y. 
Yin, and S.E. Harris, Opt. Lett. 24, 625-627 (1999); 
IEEE Journal of Selected Topics in Quantum Electronics 
5, 1502 (1999). 

[8] L.P. Yatsenko, B.W. Shore, T. Halfmann, K. Bergmann, 
and A. Vardi, Phys. Rev. A 60, R4237 (1999); T. Rickes, 
L.P. Yatsenko, S. Steuerwald, T. Halfmann, B.W. Shore, 
N.V. Vitanov, and K. Bergmann, J. Chem. Phys. 113, 
534 (2000); L.P. Yatsenko, N.V. Vitanov, B.W. Shore, 
T. Rickes, and K. Bergmann, Opt. Commun., to be pub- 
lished (2002). 

[9] R.W. Boyd, Nonlinear optics, (Academic Press, San 
Diego, 1992). 

[10] A. D. Wilson-Gordon, R. Klimovsky-Barid, and H. Fried- 
mann, Phys. Rev. A 25, 1580 (1982); L. J. Rothberg and 



N. Bloembergen, Phys. Rev. A 30, 820 (1984). 
[11] A.O. Melikyan and S.G. Saakyan, Zh. Exp.Teor. Fiz. 76, 

1530 (1979) [Sov. Phys. JETP 49, 776 (1979)]. 
[12] A.R. Karapetyan and B.V. Kryzhanovskii, Zh. Exp. Teor. 

Fiz. 99, 1103 (1991) [Sov. Phys. JETP 72, 613 (1991)]; 

B. Kryzhanovsky and B. Glushko, Phys. Rev. A 45, 4979 

(1992). 

[13] C.J. McKinstrie, G.G. Luther, and S.H. Batha, J. Opt. 

Soc. Amer. B 7, 340 (1990). 
[14] G. Cappellini and S. Trilo, J. Opt. Soc. Amer. B 8, 824 

(1991). 

[15] E.A. Korsunsky, T. Halfmann, J. Marangos, M. Fleis- 
chhauer, and K. Bergmann, unpublished. 

[16] K. Hakuta, L. Marmet, and B.P. Stoicheff, Phys. Rev. 
Lett. 66, 596 (1991); D.Kim, C.S. MuUin, and Y.R. Shen, 
Appl. Phys. B 60, S215 (1995) and references therein. 

[17] L.D. Landau and E.M. Lifshitz, Electrodynamics of Con- 
tinuous Media, (Pergamon Press, New York 1975). 

[18] P.F. Byrd and M.D. Friedman, Handbook of Elliptic In- 
tegrals for Engineers and Scientists (Springer, Berlin, 
1971); Handbook of Mathematical Functions, ed. by M. 
Abramowitz and LA. Stegun (Dover, New York, 1965). 

[19] R.B. Miles and S.E. Harris, IEEE J. Quantum Electron. 
QE-9, 470, 1973. 

[20] Similar spatial evolution of an atomic dark superposition 
has been also found for EIT-assisted frequency conversion 
in closed A- and double-A systems in the continuous- wave 
regime: E.A. Korsunsky and D.V. Kosachiov, Phys. Rev. 
A 60, 4996 (1999); J. Opt. Soc. Amer. B 17, 1405 (2000); 
D.V. Kosachiov and E.A. Korsunsky, Eur. Phys. J. D 11, 
457 (2000). 



